Recovering Unitarity of Lee Model in Bi-Orthogonal Basis 
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We study how to recover the unitarity of Lee model with the help of bi-orthogonal basis approach, 
when the physical coupling constant in renormalization exceeds its critical value, so that the Lee's 
Hamiltonian is non-Hermitian with respect to the conventional inner product. In a very natural 
fashion, our systematic approach based on bi-orthogonal basis leads to an elegant definition of inner 
product with a non-trivial metric, which can overcome all the previous problems in Lee model, 
such as non-Hermiticity of the Hamiltonian, the negative norm, the negative probability and the 
non-unitarity of the scattering matrix. 
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I. INTRODUCTION 

In 1954, T.D. Lee introduced an exactly soluble model 
to display the necessity of renormalization in quantum 
field theory (QFT) fT*!. It is delicate that with the Lee 
model, the renormalizations of the mass, wave-function 
and coupling constant were performed in a closed form 
It is worth remanding that, in the single excita- 
tion subspace, the Lee model can be reduced into the 
Fano model in atom-molecular physics [sj or the Ander- 
son model of on-site interaction free in condensed matter 
physics [J. With this recognition, the quantum manip- 
ulation for the coherent transport in quantum dot array 
has been investigated based on the Lee model [5]. 

Originally, the Lee model describes the reaction pro- 
cesses V N -\- 9, namely, a fermion V transforms to 
another fermion N by emitting bosonic 0-particle, and 
vice versa. In the single particle subspace — the V/N9 
sector, the Lee model is exactly solvable. According to 
the analytic solution, we only need to consider the mass 
renormalization of the l^-particle and coupling constant 
renormalization since the bare single-particle states of 9 
and N are both the eigenstates of the Hamiltonian, thus 
the mass renormalizations of 9 and A"-particles are not 
need, and their physical masses equal to the bare masses. 
In the V/N9 sector, the eigenstates of the Hamiltonian 
contain the physical V-stdXe and A^0-scattering states. 

In the Lee model, the renormalization about ^-particle 
results in some perplexing and interesting natures, e.g., 
the emergence of the ghost state with an negative norm. 
Actually when the physical coupling constant g, which 
obtained from the standard procedure of renormaliza- 
tion, is strong enough to exceed a critical coupling con- 
stant ^c, the bare coupling constant in the original Lee 
model is imaginary, and then the Hamiltonian becomes 
non-Hermitian with respect to the conventional inner 
product in quantum mechanics (QM) and QFT. To avoid 
this non-Hermiticity, one could restrict g to be smaller 
than (7c, but the vanishing critical value gc due to infinite 
cutoff makes g always greater than g^. In this sense, the 
non-Hermiticity is so intrinsic, that it is ineluctable. 

To overcome the ghost state problem due to the non- 



Hermiticity, a significant attempt 0, 0] is to construct 
a Hermitian, but non-local Hamiltonian corresponding 
to the modified ^-particle-Green's function of the Lee 
model by eliminating the ghost pole artificially. It was 
Kallen and Pauli [8] who first introduced an indefinite 
metric such that the norm of the physical l/-state is pos- 
itive, but the norm of the ghost state is still neg ative. 
In 1968 and 1969, Lee himself and Wick [1, [l^ used 
the same indefinite metric to discuss the unitarity of the 
^-matrix in the N9 and N99 sectors for the imaginary 
physical masses of the physical V-siaie and ghost state 
in the regime g > gc- Recently, Bender et al. [llj | intro- 
duced a different inner product by a coupling dependent 
CVT transformation to insure the positive norms of the 
all eigenstates of the non-Hermitian Hamiltonian. An 
equivalent Hamiltonian 12] for the Lee model was found 
by the similar transformation. 

In this paper, we use the bi-orthogonal basis approach 
[isl . [T^ to find a non-trivial metric and construct a new 
inner product for the Lee model in the strong coupling 
regime g > gc- In a natural and consistent way, this ap- 
proach overcome the overall problems in applications of 
the Lee model due to the non-Hermiticity of the Hamil- 
tonian, such as the negative norm and ghost state, the 
negative probability, and the non-unitarity of the scat- 
tering matrix. 

Our approach for the Lee model based on the bi- 
orthogonal basis is to use the two complete eigenstate 
sets {|i?„)} and of the Lee's Hamiltonian H and 

its conjugate . The non-trivial metric is defined by an 
operator r] (l3 |: \Dn) = r]\En), which result a new inner 
product ((/), ip) = {(j)\r]\ip). We can explicitly calculate the 
metric operator rj for both the QM Lee model (or called 
one boson mode Lee model where the 0-particle only pos- 
sesses a single mode) and the QFT model. Using this new 
inner product, we show that in the regime g > g^ the 
Hamiltonian is Hermitian, all eigenstates have positive 
norms and the scattering matrix is unitary. It is more 
important that, our obtained metric for the Lee model 
is different from that in Rcf. |lli] , and automatically in- 
sures the orthogonality of the different eigenstates. In 
this inner product space the Hermiticity of the Hamil- 



tonian implies the unitarity of the evolution. Therefore, 
with arbitrary interaction strength the QFT Lee model 
becomes acceptable in physics, and thus could be applied 
to practical systems, such as a coupled resonator array 
interacting with a two level system |5|]. 

The paper is organized as follows: In Sec. II, we re- 
consider the unitarity breaking by briefly reviewing the 
Lee model. In Sec. Ill, we introduce the concept of bi- 
orthogonal basis and use it to deal with a simple model, 
the QM Lee model, whose Hamiltonian is non-Hermitian 
with respect to the conventional inner product. We show 
find the new metric so that the Hamiltonian becomes 
Hermitian with respect to the inner product defined by 
this metric. In Sec. IV, we use the bi-orthogonal basis to 
find the new metric for the QFT Lee model when g > gc 
so that the Hamiltonian is Hermitian and the S'-matrix 
is unitary in the Hilbert space with this inner product. 
In Sec. V, the results are summarized with some com- 
ments. In Appendix, we discuss the relation between the 
Lee model and an standard model in quantum optics and 
cavity QED. 



II. BREAKING OF UNITARITY OF LEE 
MODEL 

In this section, we revisit the breaking of unitarity by 
briefly reviewing the basic properties of the Lee model 
in the conventional representation. Then we show why 
the bi-orthogonal basis is indeed needed to recovery the 
unitarity of the Lee model. 



A. Lee model 

To describe the reaction processes V ^ N + 6, the Lee 
model uses Hamiltonian H = Hq + Hi , where 



and 
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Here, V (N) is the annihilation operator of the V (N)- 
particle with bare mass mo (rriN). The operator ak (a|,) 
is the annihilation (creation) operator of the massive 9- 
particle with the dispersion relation ujk — \/k'^ + and 
the mass /z. Here, we introduce the effective mode vol- 
ume r2, the real cutoff function fk and the bare coupling 
constant go. In the infinite cutoff, namely, the function 
fk tends to unit. Here, we point that the Lee model 
neglects the processes V + 9 ^ N when the condition 
mo > rriN is satisfied. 

We focus on the V/N9 sector, i.e., the subspace with 
nv + ne = I spanned by the basis |0) and N^al |0), 





FIG. 1: (Color online) The schematic for the function h{E): 
the red curves denote the function h{E), and the solutions 
E — raM + L^k are the energy of the scattering states, (a) when 
g < Qc, the equation h{E) = possesses only one real root 
E — mv in the regime E < rriM + when h{mM -\- ii) > 0; 
(b) when g > gc, the equation h{E) = possesses two real 
roots E = mv and E — mc in the regime E < niN + uJk when 
mo < mjv + fi. 

where ny = V^V and ng = X^fc^I'^fe- ^^i^ subspace, 
we assume the eigenstates 



\E) = [cV^ + Y,<l>E{p)N^al]\Q) 



(3) 



satisfy the eigen-equation H \E) = E\E), where E de- 
note the corresponding eigenvalues. It follows from the 
eigen-equation that 



moc 
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{E - niN - ujp)(t)E{p) 
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determines the parameters c, the wave-functions (psip) 
and the eigen-energies E. Non- vanishing solutions of Eqs. 
(HI and (O lead to the secular equation h{E) = 0, where 



9V'p 



h{E) =E-mo 



1 



2iUpi} E — TTlN — t^p 



(6) 



If go is real and h{mN + fi) > (see Fig. la), there are 
two kinds of real solutions of Eq. ©, where one solution 
satisfies E < uin + ^, the others are Ek — mjq + ujk- 
Here, the condition h{mN -|- /x) > 0, i.e.. 



rriN + /i + ^ 



9lPp 



2ujpVl LOp — jl 



> mo, 



(7) 



ensures the existence of the stable F-particle |l5| whose 
eigen-energy has not imaginary part (see the appendix). 

We find physical T^-state with E = mv in the regime 
my < mN + yU as 



|V)=zy^[yt+^0^(p)7vtat]|o) 



(8) 



where the momentum representation of wave-function is 
90 fp 1 



y/2ujpn mv - mN - ojp 



(9) 
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and the normalization constant is 



(10) 



When go 0, the eigenstate |V) becomes V'^ |0). Thus, 
the dressed state |V) describes the renormahzed V- 
particle, or called the physical F-particle. The physical 
mass of the F-particle is ray, which satisfies h{'my) — 0. 
Next, we consider the scattering states 

|N0fe)± = \cu,±V^ + i>K±(.V)^^^l\ |0) , (11) 



to satisfy h{E = Ek) = for Ek = uin + LOk- While we 
assume 



gofkh-\mN + ujk ± iO+)/V^^A (12) 



the scattering features are reflected by the S-function in 
the wave-function 



0fe,±(p) =S{k-p) + 



gafpCk,: 



1 



y^2ujprt uu — ujp± ^0^ 



(13) 



The positive infinitesimal 0+ is introduced here such 
that the eigenstates ([TT]) are rewritten as the standard 
Lippmann- Schwinger scattering states 



|N0fe>± = (1 



1 



Ei,-H± iO- 



-Hi)N^ai |0) , (14) 



N^al |0) 



1 



Ek - Ho± iO+ 



Obviously, the eigenstates (|lip as well as Eq. (fTi|) de- 
scribe the scattering of the 0-particle with momentum k 
by the A'^-particle. The corresponding S-matrix element 



{N9p \mk)_ 



(15) 



is determined by the scattering phase shift = 
arctan/3fe, where 



13k 



47r(TOy — mAf — Wfc) (27r)3 
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and 



(Pp{mv - ruN - uJk) 
2u}p{ujp - ujk)imv - mN - 'jJp)'^ 



(16) 



(17) 



Here, V denotes the principal-value integral. Because 
in Eq. (fTU]) contains the divergent integral 4>v(p) 
when the cutoff tends to infinite, i.e., /p — > 1, the normal- 
ization constant Zy tends to zero so that the phase shift 
5k and the cross section both vanish. For avoiding the 
vanishing of the cross section, we introduce the physical 
coupling constant = g^Zy and assigne it to be finite. 



Breaking of Unitarity and Emerging Ghost 
State 



The normalization constant can be rewritten as 

J2 



Zy = 1 



9t 



(18) 



in terms of the critical coupling constant g^ defined as 

f2 



1 



p 



2L0p^l {my — TTlN — i^pY 



Then the relation between g^ and is 



9l 



g'gl 

9l-9^ 



(19) 
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Obviously, if g^ < g^, the normalization constant Zy 
and the square g^ of the bare coupling constant are al- 
ways positive, so that the conventional approach based on 
QFT is self-consistent and proper. However, if > g^, 
the normalization constant is negative, namely, the phys- 
ical V^-state has the negative norm. In addition, the 
square f/g of bare coupling constant becomes negative, 
so that ga = iAq (Aq G real number) is imaginary, which 
results in that the Lee model Hamiltonian is not Hermi- 
tian with respect to the conventional inner product and 
the unitarity is broken. 

As illustrated in Fig. lb, when g^ is imaginary, the 
detail analysis shows that when mg > rriN + n, there is 
no real solution satisfied E < mj\j + ^, but when mo < 
totv -I- a* there are two real solutions satisfied E < mN + fj,. 
In this paper, we consider the case mg < uin + and 
prove that ii go is imaginary, the secular equation h{E) — 
possesses another real solution E = ma (< rny) in the 
regime E < niN + H besides the real root E = my. 

To consider the physics of this solution we rewrite h(E) 
as h{E) = Zy^{E - my)F{E), where 



FiE) = 1 - 



9 



Sp{E 



my] 



{2itY J 2ujpA{E, Wp) A2 {my ,ujp)' 



(21) 



and A{E,ujp) = E — mN — tOp. The schematics of the 
functions h{E) and F{E) are shown in Fig. la (b) and 
Fig. 2a (b) for the case g > gc {g < gc), which display 
the relations dEF{E) > and F{my ) — 1 explicitly. 

It is remarkable that as i5 ^ — oo, F{E) Zy. Then 
we conclude that if g > g^ i.e., Zy < 0, there always 
exists the real root mo of equation h{E) in the interval 
{—CO, my) when mg < mN + /i, which is the energy of 
the ghost state. Here, the un-normalized ghost state is 



\G) = [V^+J2'I^G{p)N^al]\0) 



(22) 



where wave-function in the momentum representation is 
90 fp 1 



^2wpf2 mo - mN - ujp 



(23) 



F{E) 
1 - - 



(a) 




(b) 



FIG. 2: (Color online) The schematic for the function F{E): 
the red curves denote the function F{E), (a) when g < g^, 
the equation F{E) = possesses no root in the regime E < 
mjv + uJk', (b) when g > gc, the equation F{E) — possesses 
one real root E — mc in the regime E < tun + i^k- 



g <gc 



(a) 



g>gc 



(b) 



The energy spectrums of the Lee model in the V/N9 
sector are displayed in Fig. 3 schematically for the case 
g> gc and g < gc- 

To overcome the difficulty of negative norm when 
g > gc, the metric ( = (—1)"^ is introduced by Kallen 
and Pauli They found that for the case g > gc, 

in the new inner product space the norms (V| ^ |V) and 
± (N9k \ C |N6'fe)_|_ were both positive, and the orthogonal 
relations ± (N6ife| ( |V) = (V| C \'!^0k)± = were also en- 
sured. However, the norm (G|C|G) of the ghost state 
was still negative and thus the ghost state phenomenon 
still exists! The negative norm of the ghost state implies 
that the metric ^ is an indefinite metric, and the ^-matrix 
satisfies (^S^C^S = 1, which is non-unitary. Recently, Ben- 
der et al. jm] introduced the different inner product by 
the CVT symmetry to insure the positive norms of the all 
eigenstates in the sector V/N9 for g > gc- In this paper, 
we use the standard bi-orthogonal basis approach to find 
the proper inner product that automatically ensures the 
orthogonality and the positive definite of the eigenstates. 



III. BI-ORTHOGONAL BASIS USED FOR ONE 
BOSON IVIODE LEE MODEL 

In this section, we first briefiy introduce the key ideas 
of the bi-orthogonal basis with its application to the one 
boson Lee model (or called QM Lee model) where the 
0-particle o nly has one mode, which has been studied 
extensively [ill, [H, [ll| . 

A. Bi-orthogonal basis approach and Metric 
operator 

In the bi-orthogonal basis approach to the non- 
Hermitian Hamiltonian H, the eigenstates and the cor- 
responding eigenvalues |e„) and En determined by the 
eigen-equation i?|e„) = -E„ |e„) can span the whole 
Hilbert space in some sense, but arbitrary two eigenstates 



FIG. 3: (Color online) The schematic for the energy spectrum 
in the V/N6 sector: when g > gc, the ghost state emerges as 
a discrete state labelled by ma- The discrete state labelled by 
mv is the physical V^-state. Here, ruN+uik denotes the energy 
of the scattering state, mv and mc denote the energies of 
the physical V^-particle state and ghost state. 



are usually not orthogonal to each other. To have a com- 
plete orthogonal basis, we need the eigenstates \dm) of 
with corresponding eigenvalues E^- Then we have 
the two sets of the basis {|e„)} and and we can 

prove that if ^ {dm |e„) = 0, and if E^ = 
{dm Icn) is nonzero. Hereafter, we define \dn) is the eigen- 
state with the eigenvalue i?* of the Hamiltonian and 
do not assume their normalization. The two sets of the 
basis 1 1 6^)1 and {|dn)} form the so-called bi-orthogonal 
basis 

With the help of bi-orthogonal basis, the completeness 
relation is 

|en) {dn\ ^ \dn) {Cn] ^ , . 

^(d„|e„) ^(e„M„) • ^ ^ 

If we define the normalized basis by \En) = 
|e„) / ^/{(^^en} and \d„) = \dn) / \/{en M„) , we have the 
generic completeness relations 

Y,\En){Dn\=Y.\Dn){En\=l, (25) 

and 

{D.n\Em) = {En\Dm)=Smn- (26) 

According to Refs. [l3|, we can define the new inner 
product 

(a>,vl/) = ($^|vl/) = (<i>|^|vl/) , (27) 
with new metric operator rj: 

\En) -^\D„) =r^\E„)- (28) 
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The Hamiltonian H under the new metric is Hermitian, 
i.e., r^H = H^^rj or = Thus, for the 

non-Hermitian Hamiltonian that satisfies -qH = H'^t], we 
always find out the metric rj by bi-orthogonal basis, so 
that in the new inner product space the Hamiltonian be- 
comes Hermitian and the all eigenstates have the positive 
norms, and the arbitrary two eigenstates are orthogonal 
to each other. 



B. One boson mode Lee model 

In this subsection, we use the bi-orthogonal basis to 
deal with the simple model, i.e., one boson mode Lee 
model or called QM Lee model. To show the main ideas 
to solve the unitarity problem of the standard Lee model 
(the QFT Lee Model), we first consider the one boson 
mode Lee model, namely, the 0-particle only has one 
mode, which has been studied extensively [111, [H, [lB| . It 
can be proved that this model is equivalent to the "stan- 
dard model" of quantum optics, the Jaynes-Cummings 
(JC) model JJ} (see the appendix). 

The one boson mode Lee model is described by the 
Hamiltonian Hs = Hq + Hi , where 

Ho^maV'iV + mNN^N + ea''a, (29) 

and 

Hi ^ goiaV'iN + h.c.). (30) 

When go = iXo is imaginary, we use the bi-orthogonal 
basis approach to reconsider the Hermiticity of Hs- 
The eigenstates of the Hamiltonian Hg in the subspace 
spanned by the basis |0) and N'^a'' |0) are 

|*±) = M±^[i\oV^ + iE± - mo)N^a^]\0) , (31) 

and the corresponding eigenvalues are 

E±^^i.s±^u^-AXl), (32) 

where the normalization constants M± are determined 
below. Here, s = tum -I- e -I- too and u = ttin + e — mo. 
Obviously, the eigenvalues are all real only if u^— 4Aq > 0. 
Here, we focus on the case u^— 4Ag > so that all energies 
are real in this sector. Below we only consider the case 
u > 2Ao, and the other case u < — 2Ao can be discussed 
in the similar manner. 

It is clear that with respect to the conventional inner 
product (<& 1 5") the two eigenstates are not orthogonal to 
each other, i.e., (^Pip |^'±) ^ 0, and the normalization 
constants are 

A/'loid = HE± - mo), (33) 

where the subscript "old" denotes the normalization con- 
stants are defined in the conventional inner product 



space. In the dual representation the eigenstates of the 
conjugate Hamiltonian i?] = Hq — Hi are 

1 4-^) = ±N^'[-tXoV^ + {E± - mo)iVtat] |o) , (34) 

with the corresponding eigenvalues — E±. It can 
be verified that (^-^j = (^-^^l = 0, and 

(*^|*±)^0. 

To overcome the non-orthogonality difficulty, we find 
out the new metric 

|*^)(*^| + |*^)(*^| , (35) 

or its explicit form in operators: 

y/u^ - 4A^ 

-2i\a{V^Na-h.c.)], (36) 

by solving the equation = ?? |5'±). Then we define 

the new inner product (<&,*) = ($|7?|vE') = ($^| VP) = 
(<j)|v[/-D^. |vl/^) — '71^) in the single particle subspace. 
The orthogonality {^^\^ ^) = (*^| *+) = and the 
positively definite norm of the eigenstates are both en- 
sured in the new inner product space defined by metric 
rj. And the Hamiltonian is also Hermitian, i.e., rjH = 
Hh]. The new normalization condition = 
\ = 1 gives 

^fl = ^u^~4Xl{E± - mo) > 0, (37) 

which are different from the normalization constants (j33p 
determined by the conventional inner product. Then we 
also have the completeness relation 

,7(|vl/+>(vl/+| + |vl/_)(vl/_|) = 1, (38) 
(|vl/+)(vl/+| + |vi,_) (*_!),; = 1. (39) 

IV. LEE MODEL WITH g > 

In this section, we study the QFT Lee model by using 
the bi-orthogonal basis approach. The problem is con- 
sidered in the subspace with ng + ny = 1 when the bare 
coupling constant go = iAo (Aq G R) is imaginary, i.e., 
g > gc- We will prove in the appendix that the QFT Lee 
model is essentially a massive JC model where the "light 
field" would be massive. 

We first consider the ghost state and physical V^-state. 
We show that in the new inner product space, the ghost 
state and physical T^-state have positive norms. Secondly, 
we consider the scattering states of the 6'-particle. We 
prove that the scattering states also have positive norms 
in the new inner product space. 

Let me first present the new metric operator 

= bV^V + Y,WkpN^Nalap 

kp 

+ Y,{skV^Nak + h.c.), (40) 

k 
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for the Lee model as our central result. Here, the pa- 
rameters b — 2S(3 — 1 and Sp — 2'Eq(I)q{p) are defined 
by 



Wkp = dkp-2EG(t)Gik)(j)Gip)- 



(41) 
(42) 

It will be proved soon as follows that the inner product 
(<&, ^t) — ($1 1] 1^') based on this metric will solve all the 
problems in the Lee model with g > gc- (1) the ghost 
state and physical F-state have positive norms; (2) the 
scattering states of the 6'-particle has the positive norms; 
(3) the orthogonality of the eigenstates is ensured; (4) 
the unitarity of the S'-matrix for the N9 NO process 
is recovered automatically. 



A. Ghost state and physical F-state 

The energies of the ghost state and physical V-si&ie 
are located out of the continuum, i.e., tol — itin —ujk<0 
for L — V and G. In this case, the ghost state and 
physical y-state are expressed as an unified form 



\l)^e]!'[v^ 



(43) 



Here, the normalized constants Si are determined by 
the bi-orthogonal basis approach below. To find out the 
proper inner product, we consider that two dual eigen- 
states 



"l/2r 



-yt+^0^(p)Art„t]|o) 



|G^) = E'f[V^-Y,Mp)N^al]\0), (44) 
p 

of have the same eigenvalues rriL- Then we define the 
inner products (V,V) = (V^ |V) , (G,G) = (G^ |G), 
(V,G) = (V^ |G) and (G, V) = (G^ |V) for the ghost 
state and physical T^-state. The normalization conditions 
(V^ |V) = (G^ |G) ^ 1 lead to the constants S^^ in 
Eq. dH]) and 



— - 1 > 0. 

.9c 



(45) 



Obviously, with respect to the new inner products, the 
conditions (V, V) = (G, G) = 1 and (V, G) = (G, V) = 
are satisfied simultaneously. 



of Hamiltonian have the same eigenvalues Ek = tum- 
LUk- The eigenstates (|46p are rewritten as 



m 



D 



N^al |0> 



1 



Ek-Ho± iO+ 



H 



D 



.(47) 

which are the standard form of the Lippmann-Schwinger 
represntation. Obviously, the eigenstates (fT4| . ((43|) . (|44l) . 
and (|47p immediately give the orthogonal relations 



±iN9j:\mp)^ 



N0„ 



\m. 



±(Ne?|V) 



(48) 



With the help of the eigenstates of the Hamiltonian H 
and the dual Hamiltonian iJ^, we find the new metric r], 
which is satisfies rj = |^''^), so that the inner product 
is ($,^) = ($17; l^*). For the eigenstates, the relations 



N9t 



77 |V) = |V^), 77 |G> = |G^), and 77 |N0fc)± 

immediately result in the metric (j40p . Under this proper 
metric, in this subspace all states always have the positive 
norms, and the arbitrary two eigenstates are orthogonal 
to each other. And the Hamiltonian satisfies 7yi/ = H^'ri, 
which implies the Hamiltonian is Hermitian under the 
new metric. 

Under the new metric, the S'-matrix elements are de- 
fined by 



Spk^ _(N0p|77|N0fe)+ = 



|N0fc) 



(49) 



The Lippmann-Schwinger formalism results in the S'- 
matrix elements as 



Spk = Skp ~ 2niS{Ek - Ep)Tpk, 



where 



Tpk — 



gofp{Q\V\mk)+ gofpCk 



yy2ujpD, 



^2ujpfl ' 



(50) 



(51) 



T- 



defines the transfer matrix, i.e 
matrix element is obtained as Spk = 
the phase shift 6k = arctan/?fc as fk 
S-matrix is unitary under the new 
that the choice of the above matric 
representation of the S-matrix and 
observable result remains the same 
ventional matric. 



■matrix, and the S- 
exp(27(5fc), where 
1. Obviously, the 
metric. It is shown 
indeed does not the 
thus the physically 
as that in the con- 



B. Scattering states and unitarity of the S-matrix 

The energies of the scattering states are Ek — ttin + 
ujk ■ The corresponding eigenstates of Hamiltonian H are 
|N6'fe)^. In the dual space, the corresponding eigenstates 

N^fc )^ = [-Cfe,±Ut + J2 hMp)^^^] |0) ' (46) 
p 



V. CONCLUSION 

Finally, let us briefly summarize our results as follows. 
For one boson mode Lee model with an imaginary cou- 
pling constant, the Hamiltonian is non-Hermitian with 
respect to the inner product defined in the conventional 
QM, but we can use bi-orthogonal basis to find a proper 
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inner product , so that the Hamiltonian becomes Hermi- 
tian with respect to the new inner product. This inner 
product insures the orthogonahty and the positive defi- 
niteness of norm of cigenstates. For the QFT Lee model 
with g > Qc, the Hamiltonian is also non-Hcrmitian in the 
inner product defined in the conventional QFT. Based on 
the bi-orthogonal basis approach the proper inner prod- 
uct is also found. The proper inner product automati- 
cally insures the Hermiticity of the Hamiltonian, the or- 
thogonality and the positive definite norms of the eigen- 
states. Physically the ghost state is killed in our approach 
so that the unitarity of the 5-matrix is recovered auto- 
matically. 

Since the Lee model was found in 1954, there were 
many struggles to find the proper metric so that the 
positive definite norms of the eigenstates were insured. 
We have pointed out that the inner product introduced 
by Pauli et al. ^] did not insure the positive norm of 
the ghost state. Our inner product is so proper that 
the orthogonality and the positive definite norms of the 
eigenstates are ensured simultaneously. And the Hamil- 
tonian becomes Hermitian in the Hilbert space with the 
new metric, which ensures the unitary evolution of the 
quantum states. Our investigations about the Lee model 
suggest that the Hermiticity of the Hamiltonian is re- 
lated to the definition of the inner product space though 
this point is very clear in mathematics. For the non- 
Hermitian Hamiltonian in one inner product space with 
the metric, it may become Hermitian in another inner 
product space with the different metric. 

In the future work, we will used the bi-orthogonal basis 
to reconsider the unitarity of the S'-matrix in the higher 
sectors, such as NOO sector d [13, E, [H, H HH. We 
would establish the new theoretical system of QFT by 
the bi-orthogonal basis, such as the new Feynman rules, 
renormalization procedure Q and reduction formalism 
[12, [13] ■ The bi-orthogonal basis can also be used to 
reconsider the strong-coupling induced unitarity problem 
for the generalized Lee Model in quantum optics d, [H, 
[Hj [13) [2I, [1^ , atom physics [1, [lij , condensed matter 
physics [3| and the high energy physics [30, ^] . 



where |e) and \g) are the excited and ground states of a 
two level atom with energy level spacing Vt. Here, (a[.) 
denotes the annihilation (creation) operator of the pho- 
ton with dispersion relation ujk = k = \Jk'^ + + ■ 
The coupling 



G 



2eo(27r) 



rD- cos Lpk: 



of the atom to photon is determined by the angle ipk 
between the atomic dipole momentum D and the electric 
field polarization vector of photon, where eg is dielectric 
constant in vacuum. 

This model can be applied to describe the spontaneous 
emission of the atom. In the single mode limit, the model 
becomes the single mode JC model 



Hjc — ^\e) {e\ + ojf^a) a 

-hGo(a|e) (5I + /I.C.), 



(A.2) 



which is just the one boson mode Lee model or QM Lee 
model discussed in Sec. HI. 

Comparing Eq. (jA.ip with the Lee model, we find 
that the multi-mode JC model is formally the massless 
Lee model with totv = 0, /i = 0, = mg, and Gk = 
gofk/^/2ujk^. However, the Lee model has very different 
nature from that of the multi-mode JC model due to 
the non-vanishing mass /i. Thus, owing to the different 
dispersion relation, the Lee model may has the stable V- 
particle state |V), but the multi-mode JC model possess 
the unstable excited state that has the finite lifetime and 
decays to the ground state by emitting photon. 

In the following discussions, we calculate the lifetime of 
the excited state of the two level atom by the scattering 
theory. Now, we do not limit the form of the dispersion 
relation Uk — w{k), but we assign it with a lower cutoff 
LOrain for any k. In this sense, we can generally study how 
the form of dispersion relation effects on the existence of 
the bound states. 

We consider the scattering states 
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fc)± = Uk^± |e) + ^Xk,±{p)al \g) 



(A.3) 



APPENDIX: RELATIONSHIP BETWEEN LEE 
MODE AND A "STANDARD MODEL" IN 
QUANTUM OPTICS 

The JC models with one mode and multi-modes can 
be regarded as a "standard model" in quantum optics or 
cavity QED, which deals with the discrete levels inter- 
acting with some continuum [13, ■ In the appendix, 
we consider the relation between the Lee model and the 
multi-mode JC model with the Hamiltonian 

Hsh = ^ |e) (e| + y^^fca^flfc 
fc 

+ J2Gk{ak\e){g\ + h.c.), (A.l) 



corresponding to |N^^/j)^ in the Lee model. Here, 

Uk,± Gkhjf^^{u;k), 



Xk,±{p) = Skp + 



GpUk^ 



LUk — UJp «0+ ' 



(A.4) 



are determined by the eigen-equation Hsh\g,k)^ = 

(^k \g,k)_j_. We define the function hsb,±{i-j) = G^{uj)^ 
where 



G±(c^) = [w-r!-^. 



Gi 



(A.5) 



is just the retarded G+{u}) or advanced G_(u;) Green's 
function of the excited state. Then the S'-matrix elements 



Skp = -{Nep\N9k)+ are 

Skp — 5kp — 2niS{uJk — u)p)GkUk,+, 
and thus the scattering phase shift is 

If there existed the bound states 

\B) - \e)+J2xB{p)al\g), 



(A.6) 
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pole 
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(a) 


(b) 



in this system for Xsip) — GpiEg — ujp) ^, the energies 
Eb < Wmin of the bound states \B) would satisfy 



(A.9) 



Thus, we conclude that the poles of the S'-matrix (jA.6[) 
or phase shift (jA.7P determine the bound state energies. 

Using Eq. (IXO]) . we discuss the existence of the bound 
states. Because duihsb{oj) > and hsbiijJ — > — oo) — cxd, 
Eq. (IA.9P possess the real roots Eb in the regime Eb < 
t^min only if hsb{i^~iJ > 0, where uj~-^ = tj,nin - 0+. 
However, in the opposite case /is6('^,7iin) < i-*^-' 



(A.IO) 



In this 



there is no real root in the regime Eb < ijJn 
case (|A.10p . Eq. (jA.9p possess the imaginary solution 
^'quasi ^ Eb — i^y in the regime Eb > ujaiin- Though the 
solution -Equasi is not the eigenvalue of the Hamiltonian 
due to the Hermiticity of the Hamiltonian, it describes 
the unstable excited state with the energy Eb and the 
lifetime I/7. For the small decay rate, the real part Eb 
and the decay rate 7 are determined by 



and 



Reh,b,+ iEb) = 0, 

lmhsb,+ {Eb) 
d^Rehsb,±{uj)\^^ 



(A.ll) 



(A.12) 



approximately Here, Rehsb,+ {Eb) ~Rehsb,-{Eb). 

On the other hand, near the pole -Equasi the Green's 
hmction G^{uj) of the excited state is 



G+(c.) 



1 



Rehsb,±{Eb) + ilmhsb,±{Eb) ' 
Z 



UJ — Eb + i^ 



(A.13) 



where Z^^ = d^Rehsb,+ {uj)\^^^^^- It is remarkable that 
when the energy of incident photon equals to Eb^ the 
phase shift 5k goes to 7r/2 so that the cross section tends 
to infinite, which exhibits a typical resonance. 



FIG. 4: (Color online) The analytic property of the Green's 
function h~^{uj) of JC model as an Lee model with vanishing 
mass of ^-particle: the Green's function possesses the branch 

nd 
les 
nd 
jle 



□ 



@ - plane 



branch cut 



• E_ = a-ir, 



FIG. 5: (Color online) The analytic property of the Green's 
function Gjoj): because the vanishing mass /i, the excited 
state energy Q is embedded in the continuum ujk and the 
atom has the finite lifetime l/7e. 



In the Lee model, Wmin regardes the mass fj, of boson, so 
that under the condition hiruN + ^J') > there always ex- 
ists the bound state |V), i.e., the physical V-particle state 
without the imaginary part of the eigen-energy. However, 
if /i(mjv + /i) < the bound state vanishes, and the sta- 
ble IZ-particle state becomes the unstable state. Thus, we 
conclude that in the Lee model, due to the non-vanishing 
mass /i, the system contains the bound state |V) only if 
/i(r7ijv + Ai) > 0. The analytic property of the Green's 
function G+(a;) in the whole complex-w plane is shown 
in Fig. 4, which exhibits that when h(mN + fi) < 0, the 
bound state vanishes and two poles E± in the upper and 
lower planes emerge, the imaginary part 7 of the E^ is 
the decay rate of the unstable I^-particle. 

For the multi-mode JC model, the mass of the pho- 
ton is zero. We except the excited state is greater than 



9 



the ground state of the atom so that the energy of the 
excited state always locates in the continuum of the pho- 
ton energy spectrum and the excited state is unstable 
with the lifetime r ~ l/7e- The renormalized energy of 
excited-state is i7r = + <50, where the Lamb shift 



(A.14) 



in the whole complex-w plane is shown in Fig. 5, which 
exhibits that the two poles E± in the upper and lower 
planes emerge and the imaginary part 7 of the is the 
spontaneous emission rate. 



is determined by the condition 'Rehsb,+ {to) — 0. When 
Lamb shift is small, the decay rate of the excited state is 



67reo 



(A.15) 



which is just the spontaneous emission rate, which was 
given in many references [H, HI, HI] . The analytic prop- 
erty of the Green's function 



G(w) = [w - - ^ ■ 



Gl 



(A.16) 
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